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ABSTRACT 

The James Webb Space Telescope (JWST) is predicted to make great advances in the 
field of exoplanet atmospheres. Its 25 m 2 mirror means that it can reach unprece¬ 
dented levels of precision in observations of transit spectra, and can thus characterise 
the atmospheres of planets orbiting stars several hundred pc away. Its coverage of 
the infrared spectral region between 0.6 and 28 pm allows the abundances of key 
molecules to be probed during the transit of a planet in front of the host star, and 
when the same planet is eclipsed constraints can be placed on its temperature struc¬ 
ture. In this work, we explore the possibility of using low-spectral-resolution obser¬ 
vations by JWS T/NIRSpec and JWS T/MIRI-LRS together to optimise wavelength 
coverage and break degeneracies in the atmospheric retrieval problem for a range of ex¬ 
oplanets from hot Jupiters to super Earths. This approach involves stitching together 
non-simultaneous observations in different wavelength regions, rendering it necessary 
to consider the effect of time-varying instrumental and astrophysical systematics. We 
present the results of a series of retrieval feasibility tests examining the effects of instru¬ 
ment systematics and star spots on the recoverability of the true atmospheric state, 
and demonstrate that correcting for these systematics is key for successful exoplanet 
science with JWST. 

Key words: Methods: data analysis - planets and satellites: atmospheres - radiative 
transfer 


1 INTRODUCTION 


The field of exoplanet atmospheric science has expanded 
rapidly since the first transit spectroscopy measure¬ 
ment of the hot J upiter HD 209458b’s atmosphere by 
ICharbonneau et"al1 (200 2). Since then, the atmospheres of 
many hot Jupiters and some smaller planets have been 
at least partially characterised using this technique; the 
presence of clouds is sugges ted by observati ons of HD 
189733b (e.g.lPont et al.ll2Q13l j. WASP-12b (e.g. ISindl2013l: 
Stevenson et al.ll2014| j and GJ 1214b (e.g. iKreidberg et all 


2014); H 2 O absorption features have been measured for sev¬ 


eral hot Jupiters us i ng the Hubble space 1 elescope / WrU 
(e.g. [peming 2 0131 : Wakeford et al.l 120131 : iMandell et al.l 
120131 : [Frame et al.l 120141 ") : thermal emission spectra and 
phase curves have been obtained for HD 189733b and HD 
209458b (e.g. Peming et al. 200 61: Charbonneau et al.ll2008l : 


iKnutson et al. 120121 : Zellem et al l 20141 ") : and the visible re 
f lection spect rum of HD 189733b has also been measured 
([Evans et al ]|H3). 


Transit spectroscopy is a technique used to infer the 
atmospheric properties of exoplanets, when the planet tran¬ 
sits or is eclipsed by its host star. During a planetary transit, 
opacity sources in the atmosphere can make the planet ap¬ 
pear larger at particular wavelengths, allowing atmospheric 
composition to be inferred. During eclipse, or secondary 
transit, the planet’s thermal emission spectrum can be ex¬ 
tracted by comparing the in- and out-of-transit fluxes, pro¬ 
viding information about atmospheric structure. The fluc¬ 
tuations in the transit depth as a function of wavelength 
are small, of order 1/10,000 for a hot Jupiter-size planet 
and as small as 1/1000,000 for Earthlike worlds. The in¬ 
frared James Webb Space Telescope (JWST) will transform 
our ability to study exoplanet atmospheres, using transit 
spectroscopy and other methods; for example, its large col¬ 
lecting area and coronographic capability will make it an ex¬ 
tremely useful tool for direct spectroscopic imaging of hot, 
young exoplanets far from their host stars, such as the HR 
8799 system. In addition, it will greatly expand the distance 
out to which spectroscopic follow up is possible for transit¬ 
ing exoplanets, due to its large collecting area, although its 
high sensitivity limits the possibility of investigating closer 
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Figure 1. Retrieval sensitivity to temperature as a function of 
pressure and wavelength, for a hot Jupiter orbiting a sun-like 
star (top) and a hot Neptune orbiting an M dwarf (bottom). The 
contours correspond to a change in observed flux ratio for a 1 K 
change in the temperature at each level in the atmosphere. The 
division between the NIRSpec and MIRI instruments is marked as 
a thin white line. MIRI has sensitivity over a very limited pressure 
range, but has higher sensitivity than NIRSpec over that range. 
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planets as their host stars may be bright enough to saturate 
one of the instruments. Given the size of the transit depth 
fluctuations that are the measurement of interest, JWST 's 
large collecting area and high sensitivity are beneficial to 
increase the signal to noise ratio. Instrument systematics 
are also likely to play a significant role, and the low pho¬ 
ton noise should make these easier to identify and therefore 
allow us to reach the required high precision. JWST is the 
first space telescope to be designed and tested in the era of 
exoplanet spectroscopy, and slight modifications were made 
to optimise the telescope for that purpose; we also have the 
opportunity to study and find ways of correcting for system¬ 
atics prior to launch, and JWST 's intended location at the 
L2 point should minimise thermal fluctuations. 

The main difficulties in transit spectroscopy have arisen 
because available instruments generally only cover a narrow 
wavelength range, meaning that spectra taken at different 
times from the ground or from space must be stitched to¬ 
gether to obtain a coherent picture of the planet’s atmo¬ 
sphere. Broad wavelength coverage is particularly useful for 
breaking degeneracies between gas ab undances and temper- 
ature structure in secondary transit ([Barstow et al.ll2Q13al . 
Figure [TJ and gas abu n dances and cloud coverage in primary 
transit ([Barstow et al.l l2013bl h because spectral properties 
of clouds vary slowly as a function of wavelength. Prob¬ 
lems arise when stitching together non-simultaneous mea¬ 
surements of transits because the recorded baseline flux of 
the star may change between observations, due to, for ex¬ 
ample, stellar activity or variation in systematics between 
instruments/within a single instrument over time. This can 
lead to obvious mismatches between overlapping spectra 


(e.g. FORS observations of GJ 1214b bv lBean et al.l [2oTlh 
or the spurious detection of atmosp heric features (e.g . the 
GJ 1214b Ks-band observation by ICroll et~ahl l201ll . ini¬ 
tially interpreted as evidence of methane absorption in a 
low mean molecular weight atmosp here but not confirme d 
by the narrow band observations of lColon &; Gaidosll2013l ). 
JWST will eliminate a great deal of the need for stitch¬ 
ing due to the relatively broad wavelength coverage of some 
of its instruments. The broadest near-simultaneous wave¬ 
length coverage with JWST can be achieved by combining a 
near-in frared instrument with the Mid-Infrared Instrument 
(MIRI. [Reike et ahll2014l ). MIRI offers two spectroscopic ob¬ 
serving modes: the low-resolu tion prism mode (LR S), which 
covers 5—12 pmm at R~100 ([Kendrew et al.ll2014h ; and the 
medium-resolution integral field spectrometer (MRS), which 
extends from 5 to 28 qmm at R~1500-3500. While the MRS 
provides greater wavelength coverage, the LRS is more suit¬ 
able for exoplanet transit spectroscopy, as it covers the 5— 
12 qm region in a single exposure. The MRS requires three 
separate exposures to cover the full range. The LRS can in 
addition be used in slitless mode which is optimised for ob¬ 
servations of bright compact sources such as exoplanet host 
stars, and our study assumes use of this mode. 

Of the near-infrared instruments, the Near-Infrared 
Camera (NIRCam) and the Near-Infrared Imager and Slit- 
less Spectrograph (NIRISS) will both be valuable for ex¬ 
oplanet observations. However, given the importance of 
wavelength coverage for our work, we focus instead on the 
Near-Infrared Spectrograph (NIRSpec), which provides the 
largest single-shot wavelength coverage of all near-infrared 
instruments with its fixed-slit low-resolution mode. Like 
MIRI, NIRSpec offers a higher-resolution integral field spec¬ 
troscopy mode, but this cannot provide the full 0.6—5.0 
qmm in a single exposure. The low-resolution prism mode 
uses fixed slits. 

Since the Near InfraRed Spectrometer (NIRSpec) and 
MIRI cannot be used simultaneously, accurate stitching of 
spectra will be important for the success of JWST tran¬ 
sit observations. NIRSpec in prism mode, covering 0.6—5 
(im, and MIRI’s Low Resolution Spectrograph in slitless 
mode, covering 5—12 qm, provide the largest wavelength 
range possible with only two separate observations and in 
this work we focus on these instrument modes. However, our 
study of instrumental offsets and starspots will have impli¬ 
cations for combining measurements from other instruments 
too. 

In Figure [I] we show the sensitivity to temperature at 
different pressure levels in the atmosphere with NIRSpec 
and MIRI, for a hot Jupiter orbiting a sun-like star (top) 
and a hot Neptune orbiting an M dwarf (bottom). At each 
wavelength, the maximum sensitivity occurs in the region 
of the atmosphere the radiation emerges from. Sensitivity 
over a wide pressure range is required to perform spectral 
inversions and obtain the atmospheric structure, and this is 
best facilitated by broad wavelength coverage. For example, 
the MIRI instrument has little sensitivity in the deep atmo¬ 
sphere, meaning that the temperature structure could not be 
fully retrieved using MIRI alone. Whilst NIRSpec does have 
sensitivity over a wide range of pressures, including MIRI 
provides further gas absorption bands, which can help to re¬ 
solve any remaining degeneracies, and in secondary transit 
the flux ratio increases as a function of wavelength so the 
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Planet 

Mass (xlO 24 kg) 

Radius (km) 

Hot Jupiter 

1800 

75000 

Hot Neptune 

180 

30000 

GJ 1214b 

37 

15355 

‘Earth’ 

6 

6378 


Table 1. Masses and radii of the model planets in this study. 


planetary signal is stronger in the MIRI range. We therefore 
consider in this paper a combination of NIRSpec and MIRI 
observations. 

As mentioned above, NIRISS and NIRCam could also 
be used for transit spectroscopy. They both cover a similar 
spectral range to that of NIRSpec, but neither can simulta¬ 
neously observe the 0 . 6 —5 pm range, which would introduce 
further complexity when stitching together observations to 
make a complete spectrum. However, these instruments do 
have potential advantages over NIRSpec as their saturation 
limits are less strigent, so they are likely to be useful for 
observing closer, brighter systems. The bright target limit 
for NIRS pec in prism m ode is a host star J-band magnitude 
of 9.75 ([Ferruit et al.ll2014h . equivalent to GJ 1214, but 
brighter targets such as HD 189733 can be reached using 
the higher resolving power modes. In addition, the higher 
resolving power mode for MIRI is likely to be useful. The 
MIRI integral field unit extends out to 28 pm at R^3000, 
which will allow high resolution observations of, for exam¬ 
ple, the 15 pm CO 2 band, but use of this mode introduces 
further difficulty with stitching as t he observab le wavelength 
range is split into three segments (IWellsll2014l ). 

In this work, we use the NEMESIS radiative transfer 
and retriev al code (llrwin et all 20081 ) and follow the method 
adopted bv lBarstow et al.l (l2013al ) to test the accuracy with 
which combined JWST /NIRSpec and MIRI observations 
can recover the true atmospheric state of a range of planets, 
described in Section [ 2 j To explore the effect of instrumental 
and astrophysical systematic noise, we simulate the effect of 
baseline flux offsets between NIRSpec and MIRI, and also 
the effect of time varying stellar activity on temporally dis¬ 
connected observations (Section 12 . 21 ) . The method and the 
NEMESIS code are briefly described in Section [3] and our 
results are presented in Section [4j with discussion of future 
work and suggestions for appropriate strategies in Section [5] 


2 MODEL PLANETS AND SYNTHETIC 
SPECTRA 

The mod el planets we consider here are similar to those 
studied in lBarstow et al.l (|2013al ) for the Exoplanet Charac¬ 
terisation Observatory (EChO) ESA M3 mission proposal. 
The bulk properties are fisted in Table [T] So far, we focus on 
Jupiter and (mini) Neptune-sized gas giant planets, which 
are considerably simpler to model as they are assumed to be 
dominated by H 2 and He. Observation o f super Earths with 
JWST /NI RSpec has been e xplor ed by iBenneke fe Seagerl 
(120121 ) and lde Wit fe Seagerl d 20131 ). and the major difficulty 
is determining the bulk atmospheric composition from trans¬ 
mission spectra since this is degenerate with temperature 
and cloud properties. We defer further exploration of this 
to future work, although we briefly discuss the possibility of 


detecting O 3 on an Earth analogue planet in orbit around a 
nearby M dwarf. _ 

As in iBarstow et al.l (|2013al ) , we use NEMESIS to gen¬ 
erate a series of synthetic spectra based on this range of 
model planets; a hot Jupiter-size planet orbiting a Sun-like 
star (radius 695000 km), a hot Neptune/Earth analogue or¬ 
biting an M dwarf (radius 97995 km, equal to that of Prox- 
ima Centauri) and a set of spectra based on the super-Earth 
GJ 1214b. The hot Jupiter and hot Neptune don’t bear any 
direct relationship to observed planets, but the hot Jupiter 
may be thought of as an HD 209458b/HD 189733b analog, 
albeit in orbit around a slightly hotter star. These planets 
both have clear atmospheres, but we include cloud in the 
GJ 1214b and Earth analog cases: the GJ 1214b cloud is 
assumed to be made of Titan tholins and has two particle 
population s with effective radii o f 0.1 and 1 pm (for further 
details see IBarstow et a l. 2013b); t he Earth ana l og clo ud is 
a water cloud parameterised as in llrwin et"al . (120141) . We 
defer a more detailed exploration of different cloudy at¬ 
mospheres to a future publication. We use the Kurucz so¬ 
lar model spectrumQfor the Sun-like star; this is located at 
25 0 pc, significan tly furt her away than the equivalent case 
in lBarstow et al.l (|2013al ). because JWST 1 s high sensitivity 
and large mirror relative to EChO m ean that NIRSpec i n 
prism mode can be easily saturated ([Ferruit et al.l 120141 ). 
PHOENIX models of the appropriate temperature, log(g) 
and metallicity are used for the P roxim a Centauri-based M 
dwarf and GJ 1214 dHusser et al.ll2013h . The hot Neptune-M 
dwarf system is located 15 pc away, and the Earth-M dwarf 
system 10 pc away. GJ 1214 is at its correct distance of 14.6 
pc. 

For each planet, we keep the T-p profile and bulk atmo¬ 
spheric composition (N 2 -O 2 for Earth, H 2 -He for all other 
planets) constant, but we vary the trace gas abundances for 
each separate secondary transit retrieval. The atmospheres 
are assumed to be well-mixed in all cases except the Earth 
analog, which h as gas a bun dance vertical profiles as mea¬ 
sured on Earth (llrwin et al.ll20lJ ). For the primary transit 
retrieval we also vary the radius at the 10 -bar pressure level, 
as primary transit retrievals are extremely sensitive to the 
precise size of the planet. In the GJ 1214b model runs the 
optical depths of the two clouds are also allowed to vary. 
Since our aim is to explore the capability of JWST to char¬ 
acterise a GJ 1214b-like planet, the variations on the GJ 
1214b model that we use are not required to be compatible 
with the existing observed spec trum, although they are all 
similar to the best-fit case from IBarstow et all (|2013bl ). 


2.1 T-p Profiles 

The T-p profiles used in this work (with the exception of 
the Earth analog case, which uses an average Earth profile) 
are based on a single-slab atmosphere model (Equation [l]). 
We assume that the atmospheric temperature decreases adi- 
abatically through the troposphere and is either isothermal 
above that level or includes a temperature inversion. The 
model troposphere extends up to 0.1 bar, the approximate 
location of the tropopause for the solar system giant planets, 
and stars at the 1 bar level. There is very little sensitivity 
to temperature below the 1 bar level (see Figure [l]) and so 
we leave the temperature below this level constant. 
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Gas 

Source 

h 2 o 

HITEMP2010 (Rothman et al. 20101 

co 2 

CDSD-1000 (Tashkun et al. 20031 

CO 

HITRAN1995 (Rothman et al. 19951 

ch 4 

STDS fWeneer &; Chamoion 19981 


H 2 /He QBor^so^^^^rornrnhold ^89; ^or^so^^^aL ^89; JBor^sa^^^rommhold ^90; JBor^so^^^al^^OT; JBor^sow ^00^) 


Table 2. Sources of gas absorption line data. 




Wavelength («m) Wavelength (am) 


Figure 2. Example spectra for the hot Jupiter orbiting a sun-like star. Left: a single secondary transit (dashed lines indicate the 
blackbody flux ratios at 2500 K, top, and 1200 K, bottom); right: a single primary transit. Dark/light grey shading indicate lcr/2cr error 
bars. 


Ttrop = T strat - rlL^ln ( Hi) (1) 

mg \p 2 J 

Ttrop and T str at are the temperatures at the bottom and top 
of the adiabatic region, T is the dry adiabatic lapse rate, 
k is the Boltzmann constant, m is the molar weight of the 
atmosphere, g is the gravitational acceleration and pi and 
P 2 are the pressures at the bottom and top of the adiabatic 
region, in this case 1 and 0.1 bar. 

Whilst these temperature profiles are influenced only 
by very basic physics, they are sufficient to test the retriev- 
ability of temperature and composition with JWST. To in¬ 
crease the range of temperature profiles tested, we include 
some with a temperature inversion (a region of tempera¬ 
ture increase with altitude) above the tropopause. For a real 
planet, this may be caused by the presence of an ultraviolet 
or visible absorber, such as ozone in the case of the Earth 
or metal oxides in the case of hot Jupiters. In the test case 
here, it is simply included to ensure that an inversion would 
be detectable should it be present. 

Throughout this work, we refer to: the input tempera¬ 
ture profile, which is the profile used to generate the noisy 
synthetic spectra; the a priori temperature profile, which is 
the profile used by NEMESIS as the starting point for the 
retrieval; and the retrieved temperature profile, which is the 
profile that NEMESIS determined to provide the best fit to 
the noisy synthetic spectrum. If the retrieval is accurate, the 
retrieved temperature profile should closely match the input 


temperature profile and should not be influenced by the a 
priori temperature profile. 


2.2 Noise Calculation 

The noise cal c ulation in this work is similar to that of 
iBarstow et al.l (l2013al ): we assume random noise at the pho¬ 
ton limit for both instruments, using the appropriate val¬ 
ues for total instrument throughput for NIRSpec and MIRI 
in the chosen o bserva tional modes. We refer the reader to 
IBarstow et al.l (l2013aT ) for details of the noise calculation. 
The instrument properties used are as follows: for NIR¬ 
Spec, we adopt the averag e transmission and q uantum ef¬ 
ficiency properties used bv iDeming et al.l (I2009I L assuming 
a detector QE of 0.8, a telescope optics efficiency of 0.88 
and a total NIRSpec optics transmission of 0.4; for MIRI- 
L RS, we use the photo n conversion efficiencies presented 
in iKendrew et al.l (120141 ), which represents the fraction of 
photons from the source that are eventually received and 
recorded by the detector. 

In addition to the photon noise calculated for these in¬ 
struments, we also include simulated offsets due to instru¬ 
ment systematics and due to star spots. For the instrument 
systematics, we calculated the spectrally-averaged photon 
noise for each instrument, and generate a series of Gaus¬ 
sian random numbers (mean 0, standard deviation 1). We 
then add a number from this series, multiplied by half the 
NIRSpec-averaged photon noise, to all NIRSpec points, and 


1 http://kurucz.harvard.edu/stars/ 
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Temperature Profiles 



Temperature (K) 


Temperature Profiles 



Temperature (K) 


Figure 3. Temperature retrievals and error envelopes for the hot 
Jupiter secondary transit; cases shown are no baseline offset be¬ 
tween NIRSpec and MIRI (black/grey crosses, lines and shading), 
offset of order of the photon noise (blue/aqua) and offset of order 
3x the photon noise (red/pink). The input temperature profile is 
shown by the solid black line. In all cases, the a priori tempera¬ 
ture profile is an isotherm at 1500 K. The retrieved temperatures 
shown are an average over 50 retrievals with different atmospheric 
compositions, and the envelope is the standard deviation of these. 


subtract the same number multiplied by half the MIRI- 
averaged photon noise from the MIRI points. This results 
in a range of models for which the NIRSpec and MIRI spec¬ 
tra are offset from each other by a random amount on the 
order of the photon noise; for half the cases, the NIRSpec 
spectrum will be shifted above the MIRI spectrum and oth¬ 
erwise it will be below the MIRI spectrum. We also repeated 
the test using 3x photon noise, and we performed this test 
for both primary and secondary transit observations. 


2.3 Star Spots 

Star spots do not have a significant effect on secondary 
transit observations, but can alter primary transit spectra. 
This occurs because, if the star has spots on its disc, the 
brightness of the transit chord may not be the same a s 
the disc-averaged brightness. Following IPont et al.l (120131 ). 
if the transit chord is more spotty than the disc average, 
then the planet radius is underestimated because the planet 
is crossing a less bright part of the star; conversely, if the 
transit chord crosses a spot-free region, the radius is over¬ 
estimated because the occulted region is brighter than the 
average. The first case can be easier to identify, as the sig¬ 
natures of spot crossing events may be seen in the tran¬ 
sit light curve, although the wavelength ranges of MIRI and 
possibly NIRSpec are too red for these events to be seen 
above the noise since the spot contrast decreases towards 
longer wavelengths. In the second case, because no spots are 
crossed, there is not necessarily any evidence that the star 
has active regions when the transit is observed, so any poten¬ 
tial errors introduced by stellar activity may be unidentified. 
Therefore, we focus on the second case here. 

For the sun-like star, we calculate the spot contrast by 
simply assuming a star temperature of 5800 K and a spot 
temperature of 5000 K. A spot temperature 800 K cooler 


Figure 4. As Figure [3] but attempts have been made to correct 
for the baseline offsets between the NIRSpec and MIRI instru¬ 
ments. The improvement in the quality of the retrieval for an 
offset of order 3x the photon noise is significant, and a smaller 
improvement can also be seen for smaller offsets. 


than the average is comparable t o the best-fit spo t tem¬ 
perature found for HD 189733b bv IPont et al.l (|2013 ). Since 
the shortest wavelength of NIRSpec occurs longward of the 
main spectral features in solar type stars, we assume that 
the star and spots can be represented by black bodies. The 
overestimation of the squared pla net :star radius rati o is then 
caculated using the formula from IPont et ahl (120131 ): 


1/a = 


f -^5000 + (1 — /)^5800 
H58OO 


( 2 ) 


where the squared radius ratio is overestimated by a fac¬ 
tor a if the fractional coverage of spots with flux R5000 is /. 
To simulate this effect, we multiply the synthetic spectrum 
by a. The presence of star spots will have a wavelength- 
dependent effect on spectra from the NIRSpec instrument, 
since the spot contrast increases rapidly towards the visible, 
but for MIRI we treat the star spot effect as uniform across 
the wavelength range for the sun-like star. 

For the M dwarf case, we can no longer approximate the 
stellar disc and spot spectra as black bodies. This is because 
at cooler temperatures infrared absorption features, espe¬ 
cially of H 2 O, become more apparent in the spectra; since 
H 2 O absorption is increased in the cooler spots, increased 
spot coverage could potentially mimic H 2 O absorption in 
the planet’s atmosphere in the case of unocculted spots. We 
therefore use PHOENIX model spectra at 3000 K for the 
stellar disc and 2300 K (the coolest available model) for the 
spot, and calculate the spot contrast as follows: 


/ flspot + (1 f') I star / Q \ 

1/a = - - - (3) 

J- star 

At wavelengths longer than 5 microns we extrapolate 
these spectra to the closest black body curve. PHOENIX 
models at the same temperatures are also used for GJ 1214, 
with the appropriate metallicity and log(g). 

One of the biggest problems likely to arise from star 
spots is spot coverage varying in the time between NIRSpec 
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and MIRI observations of the same object. To simulate the 
effect of this, spots are applied at random to NIRSpec or 
MIRI observations, or both, or neither. This allows us to 
explore not only the effect of spots within a single obser¬ 
vation, but also how activity variation can affect stitching 
of observations. For the M dwarf stars, which we expect to 
be more active, we test a spot coverage of 10%, which is 
likely to represent the worst-case scenario. For the Sun-like 
star, the maximum spot coverage considered is 3%, as in 
general these stars are less active than M dwarf stars. How¬ 
ever, stellar activity levels for an individual star may vary 
significantly from the average for a given spectral type, so 
an activity study of interesting planet host stars before the 
JWST launch would be extremely useful. 


2.4 Offset Corrections 

As well as exploring the effect of instrumental offsets and 
star spots, we also begin to investigate how these effects 
might be corrected for. Since NIRSpec and MIRI overlap 
over a small range of wavelengths, albeit in a part of the 
MIRI spectrum where the throughput is low, it may be pos¬ 
sible to tell that an offset is present between the spectra 
measured by different instruments. If the shift is uniform, 
as assumed in the cases we test, it may also be possible to 
correct for it. 

To test this, we use NIRSpec and MIRI spectra with 
a uniform offset applied. We average the transit depth for 
each instrument in the overlapping wavelength region and 
calculate the difference. We then correct the MIRI spectrum 
relative to the NIRSpec spectrum according to the measured 
offset, and see if this improves the retrieval accuracy over the 
uncorrected case. 

We have not considered here instrumental systematics 
that result in a varying baseline across a given instrument. 
It is to be hoped that careful detector calibration and a 
detailed understanding of detector behaviour would allow 
these effects to be mitigated. 

Without monitoring the host star, a s has b een done for 
HD 189733b in the work of IPont et al.l (l2013h . it will not 
be possible to differentiate between an offset due to instru¬ 
mental effects or due to stellar activity in primary transit 
observations. Therefore, the correction applied only assumes 
a wavelength-invariant shift of the transit depth seen in one 
instrument compared to the other. This will not account for 
the non-grey effect on starspots on the NIRSpec spectrum. 


3 CODE 


We use the NEMESIS radiative transfer and retrieval code, 
which was originally developed for solar system planets and 
has since been expanded for exoplanet modelling applica- 
tion s. For deta il s of i ts use for real exo p lanet o bs ervations, 
see iLee et al. (l2012h : Barstow et al. (l2013bh: iLee et al.l 
(l2013l j: fearstow et al ] (l2014h and lLee et al ] (120131) . 

NEMESIS incorporates a fast, correlated-k based ra¬ 
diative transfer model with an optimal estimation retrieval 
scheme. Whilst optimal estimation is limited by its Gaus¬ 
sian treatment of errors, in cases where observed spectra 
contain sufficient information such schemes perform as well 
as more rigorous, but slower, algorithms such as Differential 


Evolution/M arkov-Chain Monte Carlo (iLine et al.l l2013fo . 
Correlated-k (iGoodv &; Yundfl989l : lLacis &; Oinaslll99lh is 
an approximation applied to gas absorption features within 
the radiative transfer calculation, and allows gas absorption 
(k) coefficients to be calculated offline and pre-tabulated, 
increasing the speed of the model run. K-tables can be used 
without significant loss of accuracy provided that the rel¬ 
ative strengths of absorption lines are well-correlated be¬ 
tween adjacent layers of the model atmosphere - e.g. the 
strongest line at the bottom of the atmosphere is also the 
strongest line in the next layer up. This approach has been 
extensively used to success fully model various atmospheres 
within the solar syst em fe.g. llrwin et al.ll2004lFletcherll201ll : 
iBarstow et al.ll2012h . 

We include collision-induced absorption due to H 2 and 
He within the model atmospheres, as well as absorption line 
data for H 2 O, CO2, CO and CH4. The sources for these 
data are listed in Table [2] In the Earth model we include 
data for H 2 O, CO2, CO, O 2 , O3, CH4 and N 2 O, all fro m 
the HITRAN 2008 database ([Rothman. L. S. et al.ll2009b . 


4 RESULTS 
4.1 Hot Jupiter 

The secondary transit case with JWST is the most straight¬ 
forward to deal with, because star spots are not likely to 
strongly affect these results. The only stitching problem we 
need to consider is instrumental offsets between NIRSpec 
and MIRI, so we discuss secondary transit first for each of 
the planetary types listed in Table [T] The advantage of sec¬ 
ondary transit observations is that the thermal flux emerg¬ 
ing from the planet over a range of pressures is measured, 
and so the temperature structure of the atmosphere can be 
inferred. 


4.1.1 Secondary eclipse 

The first test case is a clear-atmosphere hot Jupiter orbit¬ 
ing a sun-like star at 250 pc. The signal-to-noise ratio for 
this case is sufficient for spectral features to be detected 
(Figure [2}, so assuming no systematic instrumental offset a 
successful retrieval of temperature and H 2 O, CO 2 and CH 4 
VMRs is possible provided each species has a volume mix¬ 
i ng rat io of at least 0.5 ppmv, as found by IBarstow et al.l 
(l2013al i for ECHO (Tabled). A good retrieval of CO is more 
difficult as this molecule has fewer spectral features in the 
wavelength region considered. However, if there is a system¬ 
atic offset between the NIRSpec and MIRI instruments the 
quality of the retrieval is seriously degraded (Figure[3]). This 
is because the two instruments probe different gases at dif¬ 
ferent pressures within the atmosphere (see Figure [T] for sen¬ 
sitivity at different pressures as a function of wavelength), 
and so the offset produces artefacts in the retrieved temper¬ 
ature profile and gas abundances by increasing or decreasing 
the flux apparently emerging from the upper relative to the 
lower atmosphere, or vice versa. An offset on the order of 
the photon noise has a noticeable but fairly small effect on 
the typical accuracy of the retrieval, especially at higher 
pressures, and a larger offset of around 3x the photon noise 
seriously reduces the accuracy everywhere, again with the 
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Figure 5. 100 retrievals of gas volume mixing ratio (VMR) and 10 bar radius for the hot Jupiter in primary transit. A scale factor on 
the input temperature profile is also retrieved. Colours correspond to reduced x 2 (x 2 divided by n, the number of measurements-number 
of free model parameters). The a priori value for each gas is 10 ppmv. 


largest effect in the deep atmosphere as the sensitivity is 
lower at higher pressures. 

However, reassuringly, implementing the simple correc¬ 
tion for the shift we describe in Section \2 .41 enormously im¬ 
proves the quality of the retrieval (Figure 2]). Whilst the 
spread of retrieved temperature profiles is still larger for off¬ 
set cases, because only the relative offset between the two 
instruments is corrected for, this type of correction should 
allow us to have confidence in combined NIRSpec and M1R1 
observations. 


4-1.2 Primary transit 

We also test the accuracy of retrieved gaseous abundances 
in primary t ra nsit fo r the model planets. As discussed in 
iBarstow et al.l (l2013al i, retrieving a continuous temperature 
profile from a primary transit spectrum is not possible, since 
the deep atmosphere is not probed and therefore there is in¬ 
sufficient information to fully constrain t he T-p profile. An¬ 
other issue, mentioned in IBarstow et al.l (l2Q13bh . is that the 
planetary radius assumed at a given pressure level has a sig¬ 
nificant effect on primary transit retrievals. Since the planet 
radius is derived from the white light transit, it is dependent 
on the atmospheric properties and so must be included as 
a free parameter in a primary transit retrieval, fn this case, 
we retrieve the radius at the 10-bar pressure level in the at¬ 
mosphere, which is the base of our model atmosphere. The 
atmosphere is assumed to be opaque at all wavelengths for 


pressures greater than this (in practice, the atmosphere is 
opaque at all pressures >1 bar, see Figure [T] for sensitivity to 
different pressures in primary transit). Increasing the 10-bar 
radius therefore has two effects; firstly, it increases the total 
absorbing area of the planet, increasing the transit radius 
in a grey sense across all wavelengths; secondly, it decreases 
the gravitational acceleration at a given pressure (as this is 
proportional to 1/r 2 ) thus increasing the atmospheric scale 
height. These two effects result in significant and complex 
degeneracy between the 10-bar radius and other parameters. 
For example, increasing the H 2 O abundance in the model 
would increase the size of H 2 O absorption features in the 
spectrum, which might trade off with a smaller 10-bar ra¬ 
dius to keep the transit depth constant; however, increasing 
the H 2 O abundance will also increase the mean molecular 
weight of the atmosphere, decreasing the scale height, and 
this might then trade off with a larger 10-bar radius which 
has the effect of increasing the scale height. 

Although a continuous T-p profile retrieval is not possi¬ 
ble in primary transit, the spectral shape is not completely 
independent of temperature structure. For the hot Jupiter 
case, where the secondary transit is easily observable, we will 
have some prior information on the expected shape of the T- 
p profile. However, the terminator T-p profile is very likely 
to be cooler than the dayside profile, and won’t necessarily 
display exactly the same ve rtical shape, as observe d for the 
hot Jupiter WASP-43b bv IStevenson et ahl (|2014l ). There¬ 
fore, we also test the effect of using a generic adiabatic T-p 
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Figure 6. As Figure [5] but this time the T-p profile is fixed to the known input. 




Figure 7. The correlation matrix for a single hot Jupiter primary 
transit retrieval (temperature is not retrieved in this case). There 
are significant negative correlations between the 10-bar radius 
offset and the gas abundances. 


profile of an appropriate effective temperature, and includ¬ 
ing a scaling factor on this profile in the retrieval. This allows 
some constraint to be placed on the stratospheric temper¬ 
ature, and should ensure an accurate retrieval of gaseous 
abundances even when the T-p profile is unknown. 

In Figure [5] we show the retrieved versus input results 
for the case where the temperature is scaled; in Figure [6] 
we show the same results for the case where input temper¬ 
ature profile is used. The scatter on the results in Figure [5] 
is much larger than seen in Figure [6j as the addition of 
the temperature parameter increases retrieval degeneracy. 


To resolve this difficulty for a hot enough t a rget, obtain¬ 
ing a full phase curve as in IStevenson et al.l (120141 ) would 
provide more direct information about the terminator tem¬ 
perature. It is clear from Figure [6] that, even where the input 
temperature is used, there is some degeneracy between the 
retrieved gas abundances and the 10-bar radius. There is sig¬ 
nificant negative correlation between 10-bar radius and the 
gas abundances, resulting in a systematic underestimation 
of the 10-bar radius and overestimation of the gas abun¬ 
dances. We show the correlation matrix from the retrieval 
in Figure [71 this is the normalised error covariance matrix 
for the atmospheric state vector, showing correlated errors 
between different variables. Significant correlations, as seen 
here, mean that the retrieval is degenerate. This can be rec¬ 
tified by observing multiple transits to reduce the noise on 
the spectrum - observing 10 transits provides tighter con¬ 
straints and reduces the degeneracy, removing these system¬ 
atic under- and over-estimations. The retrieval accuracies for 
a single transit are listed in Table [H 

As for the secondary transit case, instrumental system- 
atics may need to be accounted for, and for primary transit 
there is also the possibility of variable stellar activity caus¬ 
ing problems. We test this by adding an unocculted star 
spot offset, as described in Section [2.21 at random to some 
NIRSpec and some MIRI spectra within the test set of 500 
models. This means that some combined spectra are spot- 
free, some have the same spot offset in both NIRSpec and 
MIRI sections, and some spectra have a spot offset in either 
the NIRSpec or MIRI section but not in both. 
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Figure 8. As Figure [2] but for single secondary/primary transits of the hot Neptune orbiting an M dwarf. The input temperature profile 
is shown in Figure [9] Dashed lines on the left indicate the blackbody flux ratios at 1500 K, top, and 750 K, bottom. We also show the 
positions of the main absorption bands on the right-hand plot as they are easy to distinguish for this case. 
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Figure 9. 100 eclipse retrievals of temperature and gas abundances for a hot Neptune orbiting an M dwarf at 15 pc distance. Colours 
correspond to % 2 /n. As for the hot Jupiter, temperature and H 2 O, CO 2 and CH 4 VMRs are retrieved accurately in most cases, apart 
from one retrieval with a high reduced x 2 which has low abundances of all gases so fewer detectable features. 


For a sun-like star, we find that even for a reasonably 
high spot coverage (3%) the quality of the retrieval is not 
significantly decreased. This is because the peak in star/spot 
contrast occurs at wavelengths shorter than JWST’s range, 
so especially for the MIRI instrument spots have very little 
effect. Baseline offsets between instruments can also largely 
be accounted for. 

For a cloud-free, well-mixed atmosphere, the gas abun¬ 


dance information retrieved from primary transit is very 
similar to that retrieved from secondary transit. However, 
primary transit measurements are still necessary, despite the 
added complications from stellar activity, because a) the at¬ 
mosphere at the terminator may be different from the day- 
side, and b) if the atmosphere is cloudy or not well-mixed 
then the information content of the two geometries is very 
different. Generally, primary transit spectra probe higher 
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Gas 

Hot Jupiter 

Hot Neptune 

h 2 o 

3x (>0.1 ppmv) 

1.5x (>1 ppmv) 

C0 2 

2x (>0.5 ppmv) 

3x (>10 ppmv) 

CO 

10x (>10 ppmv) 

5x (>10 ppmv) 

ch 4 

3x (>0.1 ppmv) 

2x (>0.1 ppmv) 


Table 3. Achievable precision in gas abundance retrievals from a 
single observation for the hot Jupiter and hot Neptune in eclipse. 
These precisions are only valid for certain concentrations, which 
ranges are given in parentheses for each case. As the gas abun¬ 
dances vary over several orders of magnitude, the precisions are 
given as multiples of the abundance; for example, the H 2 O VMR 
can be retrieved to within a factor 3 of the true value for a hot 
Jupiter, provided there is more than 0.1 ppmv present. The val¬ 
ues quoted take into account the errors due to retrieval degen¬ 
eracy, which is the dominant source of error, by calculating the 
precision based on the 2cf deviation from the true value over all 
retrievals. These are calculated assuming no systematic offsets. 
The hot Neptune case has a higher signal to noise as the system 
is significantly closer and the planet:star radius ratio is larger. 


Gas 

Hot Jupiter 

Hot Neptune 

GJ 1214b 

h 2 o 

3x (>10 ppmv) 

1.5x (>0.1 ppmv) 

5x (>50 ppmv) 

C0 2 

2x (>1 ppmv) 

1.5x (>0.1 ppmv) 

5x (>50 ppmv) 

CO 

Not detectable 

1.5x (>10 ppmv) 

N/A 

ch 4 

3x (>5 ppmv) 

1.5x (>0.1 ppmv) 

5x (>50 ppmv) 


Table 4. Information available from a single observation for the 
hot Jupiter, hot Neptune and GJ 1214b in transmission. Perfect 
knowledge of the temperature profile (and cloud top altitude in 
the case of GJ 1214b) has been assumed. Otherwise everything is 
as Table [3] 


regions of the atmosphere due to the longer path length in 
limb geometry, so for a case where photochemistry alters 
the composition of the atmosphere at higher altitudes there 
might be a wealth of extra information to be gained, and 
of course primary transit spectra are excellent indicators of 
the presence of clouds. 

4.2 Hot Neptune 

4-2.1 Secondary eclipse 

As in lBarstow et al.l (l2013al ) we consider a hot Neptune in a 
close orbit around a cool M dwarf. In this work, the system 
is placed at 15 p c, almost thre e times further away than 
the equivalent in iBarstow et al .1 (|2013ah , to avoid the risk 
of saturation. Spectra are shown in Figure [8] With a single 
secondary transit observed with each instrument a good con¬ 
straint can be obtained on H 2 O, CO 2 and CH 4 abundances, 
as well as the tropospheric temperature profile (Figure [9]). 
Again, a retrieval of CO is more difficult than retrievals of 
other species due to the smaller number of CO lines. 

We also find that introducing an offset between NIR- 
Spec and MIRI has a similar effect for the hot Neptune as 
we find with the hot Jupiter. Assuming there is sufficient 
overlap between the two instruments, applying a shift to 
correct for the offset means that the quality of the retrieval 
is not noticeably degraded. For secondary transit, with suffi¬ 


ciently high SNR, systematic effects such as this should not 
result in problems. 


4-2.2 Primary transit 

The effect of star spots is most critical for planets orbiting 
M dwarfs. These cooler stars generally have higher activity, 
and because of their lower temperatures their spectra con¬ 
tain some molecular absorption features, particularly due 
to water vapour. These spectral features will have different 
shapes within the cooler spot regions compared with the 
rest of the disc, and therefore unocculted spots could not 
only lead to an overestimation of the radius ratio at shorter 
wavelengths, but also inaccurate retrieved values of water 
vapour abundance. 

We present in Figure [TO] the effect of 10% spot coverage 
on the retrieval of H 2 O and CO abundances from a single 
primary transit of the hot Neptune. The retrievals for the 
spot-free, no-offset case (row 1) show that H 2 O can be re¬ 
trieved to very high precision, with CO well retrieved as 
well for abundances down to 10 ppmv. However, if instead 
the star has 10% out-of-transit-chord spot coverage for some 
of the observations, it is clear that in many cases the H 2 O 
abundance is overestimated by up to an order of magnitude 
(row 2). This is due to the differing water vapour signatures 
in the stellar and spot spectra mimicking increased water 
vapour absorption in the planet’s atmosphere. Correcting 
for the offset between NIRSpec and MIRI (row 3) doesn’t 
solve this problem; however, it does improve the CO retrieval 
from the uncorrected spot case. We can see from this that 
for CO the offset between the two instruments is the most 
important factor, whereas this is not the case for H 2 O. This 
result shows that stellar activity may be the limiting fac¬ 
tor for accurate retrievals of JWST observations for planets 
around active stars. 

The situation presented above is a worst-case scenario, 
since it may be possible to estimate the effect of star spots by 
observing spot crossing events duri ng transit s and e xploiting 
JWST ’s wavelength coverage. iFraine et al ] (l2014h use evi¬ 
dence of spot crossings in simultaneous Kepler and Spitzer 
lightcurves of the warm Neptune HAT-P-11 b to place con¬ 
straints on the temperature contrast of the star spots. Spot 
crossings are clearly observed in the Kepler lightcurves but 
not at all in the 4.5 qm Spitzer data and only indirectly in 
the 3.6 qm data, so a lower limit on the spot temperature 
can be obtained. In this case, even the coolest spots do not 
present a strong enough contrast to mimic the H 2 O absorp¬ 
tion seen in the planet’s atmosphere. It will not be possible 
to obtain simultaneous observations of spot crossings over 
JWST ’s full wavelength range, but observation of several 
transits with both NIRSpec and MIRI could provide a sta¬ 
tistical estimate of the average spot temperature, assuming 
that the planet’s orbit crosses the spot latitudes on the star. 
If this is not the case, then scenarios like the one described 
above could still be problematic, although as close-in plan¬ 
ets generally have misaligned orbits some spot crossings are 
likely. 



















Transit spectroscopy with JWST: Systematics 11 


H20 


-o 

CD 

Q) 

cc 



Input (V.M.R) 
H20 


H20 


z/n 


CO 


3.0 


2.5 


2.0 


10 

10 

10 


3 10 
> 


1 o’ 2 1-5 


I 


10 

10 

10 



//n 



120 

100 

80 

60 

40 

20 




cc 

CD 

> 

Q) 

CD 

DC 



//n 


40 

30 

20 

10 




-o 

CD 

> 

© 

'i— 

CD 

DC 




CO 


CO 




Figure 10. The effect of star spots on a hot Neptune primary transit observation. Top row: spot-free case. Middle row: 10% spot coverage 
in some observations. Bottom row: 10% spot coverage but shift between NIRSpec and MIRI corrected. 


4.3 GJ 1214b 

4-3.1 Secondary eclipse 

No observation of GJ 1214b in secondary transit has yet 
been possible, as current instruments are not sensitive 
enough to detect it since the planet is relatively cool. We 
model here what a secondary eclipse of GJ 1214b might 
be expected to look like (Figure fill). Whilst the secondary 
eclipse is certainly detectable with JWST and the broad 
shape of the spectrum is discernible, it is likely that the pres¬ 
ence of clouds on the planet coupled with a low SNR will 
confuse interpretation of the secondary transit spectrum. 
This means that retrievals of gaseous abundances and a de¬ 
tailed retrieved temperature profile from secondary transit 


are unlikely. In addition, very little flux is observable coming 
from the planet in the NIRSpec range, so just using MIRI 
is likely to be adequate and therefore systematic offsets be¬ 
tween the instruments are not a concern here. 

The most useful piece of information we could gain from 
a secondary transit of GJ 1214b is a constraint on the strato¬ 
spheric temperature, since the MIRI weighting functions 
peak in the stratsophere. Allowing for likely variation be¬ 
tween the dayside and the limb, this could be used as prior 
information for primary transit retrievals. It could help to 
break degeneracies between the atmospheric temperature, 
mean molecular weight and 10-bar radius parameter, all of 
which affect the scale height of the atmosphere. 

For 15 eclipses of GJ 1214b, a secondary transit spec- 
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Figure 11. As for Figure [2] this is a secondary transit spectrum 
for 15 eclipses of GJ 1214b. No spectral features are visible above 
the noise level. Dashed lines indicate the blackbody flux ratios at 
800 K, top, and 450 K, bottom. 



Temperature (K) 

Figure 12. An average of six temperature retrievals from 15 
eclipses of GJ 1214b, starting from six different isothermal priors 
between 400 and 650 K. All retrievals converge on the correct 
temperature at around 0.01 bar, within the grey shaded area. 
The purple line is the average, the pink shading indicates the 
spread of different retrieved values and the solid black line is the 
input temperature profile. 

trum can be obtained with sufficiently high SNR to constrain 
the stratospheric temperature at 0.01 bar by fitting the near- 
blackbody radiance (Figure fl2l) . Whilst there is clearly not 
enough information in the spectrum to retrieve the full pro¬ 
file, due to the lack of spectral features, in primary transit 
regions of the atmosphere below this are not probed anyway, 
so the detailed temperature structure at lower altitudes is 
not relevant to primary transit retrievals. 

4-3.2 Primary transit 

GJ 1214b has already been extensively observed in primary 
transit, but over the limited wavelength range available in 
space and from the ground (less than 5 pm) since Spitzer 
lost its cryogenic cooling capability. GJ 1214b has posed a 
challenge for transit spectroscopy due to its extremely flat 


spectrum, which must b e due to the presence o f cloud or haze 
high in its atmosphere ([Kreidberg et al.ll2014h . Because the 
spectrum as currently observed is featureless, no conclusions 
can be drawn about the composition of the atmosphere. Re¬ 
ducing the noise on the spectrum and increasing the wave¬ 
length coverage should help to solve this problem. 

We present retrievals (temperature and cloud top pres¬ 
sure are fixed at the correct input values) for a single tran¬ 
sit of GJ 1214 b, assum ing cloud properties as modelled in 
iBarstow et al.l (l2013bh . Values for H 2 O, CO 2 and CFG can 
be accurately retrieved to within a factor of 3 if the abun¬ 
dance of any of these gases is above 10 ppmv (Figure [13]), 
because only if they are sufficiently abundant can their spec¬ 
tral features be seen above the cloud deck (1 mbar in this 
case). The cloud optical depths and radius offset are not well 
retrieved, due to high correlations between variables (Fig¬ 
ure 0 . For a detailed discussion of degeneracies between 
cloud top pres sure and tem pe rature in primary transit ob¬ 
servations, see IBarstow et all (l2013bh . 

As for the hot Neptune case, star spots introduce fur¬ 
ther difficulties with gas retrievals for GJ 1214b, although 
the effect is less noticeable due to the relatively poor quality 
of the retrieval for the spot free case. If the spot coverage is 
up to 10%, water vapour abundance can be overestimated 
by as much as an order of magnitude (Figure [15]), and un¬ 
corrected offsets between NIRSpec and MIRI also increase 
the scatter on all retrieved gas abundances. 

Whilst JWST represents an enormous improvement on 
current instrumentation for a planet like GJ 1214b, atmo¬ 
spheric retrievals are still not straightforward due to the 
cloudiness of the atmosphere. The effect of the clouds is to 
obscure gaseous absorption features, thus reducing the in¬ 
formation content of the spectrum and making it more diffi¬ 
cult to break degeneracies in retrieval. This can be improved 
upon further by combining the results of multiple observa¬ 
tions to increase the SNR of any atmospheric features that 
are visible, but ultimately the extent of the cloud coverage 
is the limiting factor for transit spectroscopy of this kind of 
planet. 

Despite these degeneracies, some constraint on the 
cloud can be obtained with JWST. The current data from 
HST /WFC3 ([Kreidberg et al.ll2014h allow an upper limit to 
be placed on the cloud top pressure, but a range of cloudy 
atmosphere models with the cloud at lower pressures are 
still permitted. In Figure [16] it is clear that over the wave¬ 
length range probed by JWST/NIRSpec+MIRI the cloud 
top pressure degeneracy starts to be broken. In this figure, 
we show noisy synthetic spectra based on two different at¬ 
mospheric models, one with a cloud top at 0.01 mbar and 
one with a cloud top at 0.1 mbar. Both models produce an 
equally good fit to the existing data, but have rather differ¬ 
ent properties. To compensate for the higher cloud top, the 
0.01 mbar model has a smaller cloud optical depth, which 
results in the extinction at longer wavelengths being smaller 
than for the 0.1 mbar model which produces a much flatter 
spectrum. In Figure [TT1 we demonstrate sensitivity to cloud 
top pressure by comparing retrievals of H 2 O abundance from 
15 transits of GJ 1214b, for the case where the input cloud 
top pressure is used in the retrieval versus the case where a 
higher cloud top pressure is used. The spectral fit is consis¬ 
tently worse if the wrong cloud top pressure is used, again 
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Figure 13. 100 retrievals for a single transit of GJ 1214b observed with JWST. Gases are only well retrieved if the abundances are 
above a 50 ppmv threshold, when the features can be distinguished above the broad signature of the cloud. Gas priors are 100 ppmv. 


demonstrating that we can use JWST observations to place 
some constraint on cloud top pressure for GJ 1214b. 


4.4 Earth 

We perform a final test in primary transit for an Earth ana¬ 
log planet in orbit around an M dwarf at 10 pc distance. 
This is a very difficult target to observe even with the sensi¬ 
tivity of a telescope like JWST, but if a close enough system 
should be discovered it is necessary to think about what we 
might be able to see. 

Ultimately, we would like to be able to detect gases 
that might indicate an Earthlike environment, or even life. 
Ozone is a key infrared absorber that is chemically linked 
to diatomic oxygen, so it would be extremely interesting 
if we could detect this in the atmosphere of a temperate, 
Earth-size planet. We simulate a spectrum for 30 transits of 
an Earth in orbit around an M dwarf, which would span of 
order one year as each orbit would take appro ximately 10 
days. We use the Earth atmosphere model from llrwin et al.l 
(120140 . Although the resultant primary transit spectrum is 
clearly very noisy indeed, it is possible to pick out the 4.3- 
|um CO 2 band and the 9.6 -\im O 3 band; these can be seen 
more clearly in the binned spectrum (Figure [151) . The spec¬ 
trum is too noisy for an accurate retrieval of the gas abun¬ 
dances, but the presence of ozone can be detected using a 
retrieval technique; even when the ozone prior is set to very 
low values, the retrieved value is similar to the true abun¬ 


dance. No offsets or star spot corrections were applied to 
this spectrum, but even in the presence of these systemat¬ 
ics we would expect this feature to remain detectable as a) 
the feature sits in the middle of the MIRI wavelength range, 
so should be relatively unaffected by instrumental baseline 
offsets and b) the effect of star spots at these wavelengths 
is smaller and relatively grey. We would also not expect to 
see an O3 feature within a stellar spectrum; H2O features 
within the stellar and spot spectra are the main cause of the 
problems retrieving H 2 O in the hot Neptune case presented 
here. 

Of course, the photochemical environment of an Earth¬ 
like planet in a close orbit around an M 5 star is very di fferent 
from the environment around the sun. iGrenfell et al.l ( 20131 ) 
find that O 3 would still be produced around cooler stars, but 
the mechanism is likely to be different and this may result 
in differences in the vertical profile, ft could therefore con¬ 
ceivably be the case that O 3 is produced, but not detected 
because it is deeper in the atmosphere and the signature is 
obscured by clouds. 


5 DISCUSSION 

JWST will clearly be a powerful observatory for transit 
spectroscopy, and undoubtedly for other methods of exo¬ 
planet characterisation not discussed here. However, there 
is a significant amount of work to be done in order to guar¬ 
antee that we can fully exploit its potential. JWST’s lifetime 
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Figure 14. The correlation matrix for a single GJ 1214b primary 
transit retrieval (temperature and cloud top pressure are not re¬ 
trieved in this case). There are significant negative correlations 
between the 10-bar radius offset and the gas abundances. 


is limited by the amount of propellant it can carry, and so it 
is especially crucial that multiple observations of small, low 
SNR targets are planned well in advance. We here identify 
several areas for necessary future work. 


5.1 Detector systematics 

We have highlighted the issues that can arise if there are 
baseline offsets between instruments for which we are stitch¬ 
ing together observations, and discussed ways in which these 
might be corrected. The effect of these offsets before correc¬ 
tion makes it clear that detector systematics are likely to be 
important in transit spectroscopy with JWST, and a base¬ 
line offset is only one simple example. Infrared detectors 
such as those in the Spitzer 1RAC ins trument can display 
intrapixel gain v ariability (seen bv e.g. iKnutson et al.ll2008l 
for HD 209458b, lAnderson et aljfeoill for WASP-17b) and 
other more complex systematic effects that require more so¬ 
phisticated techniques to correct for than those discussed in 
this wo rk. Methods such as independent com ponent analy - 
sis (e.g. Waldmann 20l3) ? Gaussian Processes (Gibson 2014) 
and pixel mapping have all been used to correct for these 
more complex effects when reducing the observed data. How¬ 
ever, for JWST we have the advantage of studying the de¬ 
tector systematic behaviour prior to launch with a view to 
using them for exoplanet spectroscopy - a luxury that was 
unavailable for HST and Spitzer. This, combined with the 
continuing development of data reduction techniques, should 
ensure that detector systematics are less of an obstacle for 
JWST. 


5.2 Starspot monitoring and activity 
characterisation 

It is clear that, particularly for cooler stars, the effects of 
stellar activity and especially unocculted spots must be ac¬ 
counted for in transit spectra. At short wavelengths, unoc¬ 
culted spots can mimic a larger planet:star radius ratio and 
also a larger water vapour abundance. This could poten¬ 
tially result in false positive detections of cloud and water 
vapour, as well as inaccurate retrievals of other quantities. 
Unfortunately, it is difficult to blindly separate the effects 



H20 



Input (V.M.R) 


x 2 /n 

.1 


ay 

//n 

4.0 

3.5 
3.0 

2.5 

2.0 

15 l 


Figure 15. The effect of 10% out-of-transit star spot coverage on 
the retrieval of water vapour for a GJ 1214b-like planet. As for 
the hot Neptune, the presence of unocculted starspots can lead 
to an overestimation of water vapour that is difficult to correct 
for. 


of stellar activity from genuine atmospheric features within 
transit spectra. 

Attempts have been made to overc ome this difficulty 
for recent observations of HD 189733b. IPont et ahl £ 20131 ) 
correct for the effects of starspots by monitoring the flux of 
the host star HD 189733 over a long period of time using the 
Automated Patrol Telescope at Siding Spring Observatory, 
which provides an indication of the varying spot coverage. 
This time series is used to estimate the likely spot coverage 
of the star at the time of various observations by HST and 
Spitzer , and the resulting radius overestimation factor due 
to unocculted spots is calculated. Occulted spots are also ac¬ 
counted for. Whilst this method still relies on estimating the 
(non-zero) spot coverage when the stellar flux is at its peak, 
it represents a vast improvement on methods that make no 
attempt to account for star spots. It also helps when stitch¬ 
ing together measurements from multiple instruments as it 
implicitly corrects for baseline offsets due to stellar activity. 

Of course, it is not possible to monitor the activity of 
every potential target star for JWST, but for the most in¬ 
teresting targets around active stars it is an approach that 
must be considered. It relies on the availability of small, 
ground-based telescopes such as APT. 
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Figure 16. Primary transit spectra based on two different atmospheric models for GJ 1214b (15 transits). Both are compatible with 
current data, but the y diverge at wavele ngths that will be sam ple d for the first time with JWST. Existing data are as shown in 
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Figure 17. Retrieval of H 2 O abundance from 15 transits of GJ 1214b if the cloud top pressure used in the retrieval is the same as the 
input (1 mbar, left) and if it is higher (0.1 mbar, right). The x 2 / n values where the incorrect cloud top pressure is used are higher, 
indicating that the fit is not as good and demonstrating sensitivity to cloud top pressure. 


5.3 Target searches 

In order to fully exploit JWST's sensitivity, we require more 
super-Earth targets orbiting relatively nearby M dwarf stars. 
Over the next few years, ground-based and space-based sur¬ 
veys should provide more of these interesti ng targets. The 
Next Generation Transit Survey (NGTS) ([Chazelas et al.l 
12012 I) is particularly suited to this because it operates be¬ 
tween 600 and 900 nm, a wavelength range closer to the 
flux peak of cooler stars. Ideally, JWST would undertake 
an atmospheric survey of a range of super-Earth planets 
similar to GJ 1214b, to investigate the typical range of at¬ 
mospheric types and conditions, but at the moment there 
are very few for which follow up is possible. JWST' s high 
sensitivity may make follow up easier for more distant tar¬ 
gets, but if the majority of these small, super-Earth worlds 
are as cloudy as GJ 1214b, atmospheric signals for planets 
around more distant stars are still likely to be prohibitively 
small. Therefore, an important goal for the few years lead¬ 
ing up to the launch of JWST is the detection of other 
super-Earth planets at comparable distances to that of GJ 
1214b. The MEarth survey that discovered GJ 12141 3 and 
its counterpart, MEarth South ([irwin. J. et al lllill) , are 
likely to contribute to JWST' s target list of small planets, as 


will t he K2 mission (|Demorv et al.l 20 13: Vanderburg et al.l 
120141 ) and SPECULOOS survey JGillon et al .l 20131). T he ad- 
dition of the HAT-South network ([Bakos et al.ll2009h to the 
radial velocity search effort should also yield more planets 
with constrained masses. 

Other surve ys such as th e Transiting Exoplanet Survey 
Satellite ( TE67?. [Rlckerll2o3 ) . due to launch in 2017 to look 
at nearby bright stars, and contin uing efforts su ch as the 

will also 


WASP survey for hot Jupiters 
provide interesting targets for JWST. 


5.4 Cloudy atmosphere models 

Clouds clearly play a very important role in the atmosphere 
of GJ 1214b, and there is also evidence of their presence in 
hot Jupiter atmospheres. Until recently, the majority of at¬ 
mospheric models have been cloud-free, but as the evidence 
for exoplanet clouds increases this is no longer a reason¬ 
able approach. JWST transit observations should aid dis¬ 
crimination between cloudy and clear atmospheres, and for 
favourable targets may allow constraints to be placed on 
the composi tion, structure and pa rticle size of any clouds 
present (e.g. lWakeford &; Sindl2014h . We defer more detailed 
explanation of this issue to a future paper, which will focus 
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Figure 18. As Figure [2] but for 30 primary transits of an Earth orbiting an M dwarf. The spectrum on the right is binned up by a factor 
5 to make the ozone band at 9.6 microns more obvious. The CO 2 band at 4.3 microns is also clearly visible. 


explicitly on determing cloud properties using JWST trans¬ 
mission spectra. 


5.5 Other instruments and observing scenarios 

In this paper, we have considered the instruments and modes 
which allow us to maximise the wavelength coverage ob¬ 
tainable with JWST for the smallest number of separate 
transit observations. Other observational scenarios are pos¬ 
sible, and may indeed be preferable for secondary transit or 
for less active stars. Two other near-infrared instruments, 
NIRISS and NIRCam, will also have the capability to do 
transit spectroscopy, although neither can capture the 0.6— 
5 qm wavelength region simultaneously. The NIRISS slit less 
spectroscopy mode covers 0.7—2.5 qm and is specifically de¬ 
signed for transit spectroscopy, so this will be a useful tool 
for capturing that specific wavelength range. NIRCam offers 
the opportunity for slit less spectroscopy between 2.5 and 5 
qm, at R 2000. 

NIRSpec and MIRI can both be used at higher resolving 
power than considered here, and whilst this also limits the 
instantaneous spectral range achievable, these modes will no 
doubt be useful. For example, higher resolving power can 
resolve the shapes of bands in secondary transit sufficiently 
to provide constraints on vertical variation of trace species, 
or probe narrow line cores above the cloud level in primary 
transit that would otherwise be missed. The MIRI integral 
field unit’s greater wavelength range provides the potential 
to place stronger constraints on gases such as CO 2 that have 
significant bands longwards of 12 qm, and may also be useful 
for secondary transits of cooler planets where their flux may 
be expected to peak. The bright target limit for NIRSpec is 
also less stringent when using higher resolution modes, as the 
photons are spread out more across the detector, avoiding 
problems of saturation. In future work we aim to consider 
some of these observations; however, we stress that stitching 
together more segments to make a coherent spectrum will 
require a better understanding of systematics. 


6 CONCLUSION 

JWST will be a powerful tool for exoplanet transit spec¬ 
troscopy of a wide range of planets, due to its broad wave¬ 
length coverage and high sensitivity. This is the first space 
telescope built in the exoplanet spectroscopy era, and so for 
the first time we have the opportunity to test and charac¬ 
terise detectors prior to launch with a view to using them 
for transit spectroscopy. This should allow us to achieve the 
required precision of 10—100 ppm in transit depth. 

JWST will allow us to characterise the temperature 
structure and atmospheric composition of hot Jupiters orbit¬ 
ing bright stars, with gas abundances retrieved down to 0.5 
ppmv even for systems at large distances, ensuring a large 
range of accessible targets. It will also improve our under¬ 
standing of cloudy super-Earths such as GJ 1214b by provid¬ 
ing more wavelength coverage to allow distinction between 
different cloud models, and would enable us to determine 
some properties of Earth-like planets around cool stars pro¬ 
vided they are close enough. However, accurate retrievals 
for many of these objects will depend on either observing 
relatively inactiv e stars, or using stellar monitoring as in 
IPont et ahl J2013l i to correct for the presence of starspots. 
For high SNR targets such as the hot Neptune-M dwarf 
system considered here, stellar activity is likely to be the 
limiting factor on how accurately we can constrain the at¬ 
mosphere, as a 10% star spot coverage can lead to errors of 
up to an order of magnitude in the H 2 O abundance retrieval. 
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